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I.  INTRODUCTION 

This  report  summarizes  the  accomplishments  on  U.  S.  Army  Research 
Office  Grant  No.  DAAG29-76-C-0067  from  15  October  1975  to  14  October 
1978.  The  purpose  of  this  three  year  effort  was  to  develop  computa- 
tional techniques  applicable  to  the  analysis  and  design  of  electrically 
small  antennas,  and  to  use  these  techniques  to  study  small  antennas. 

Three  separate  topics  were  considered,  each  of  approximately  one 
year  duration.  Briefly,  they  are: 

1)  1st  year  - wire  antennas  in  the  presence  of  material  bodies 

2)  2nd  year  - closely  spaced  thin  wires 

3)  3rd  year  - small  antenna  location  synthesis. 

Each  of  these  topics  resulted  in  a journal  publication  and  this 
report  will  summarize  the  work  and  indicate  its  relevance  to  small 
antennas.  The  full  journal  publications  are  included  in  the  Appendixes 
for  those  desiring  additional  details.  The  summaries  in  this  report 
also  list  oral  publications,  personnel  involved,  and  degrees  granted. 

II.  FIRST  YEAR  - WIRE  ANTENNAS  IN  THE  PRESENCE  OF  MATERIAL  BODIES 


Journal  Publication  E.  H.  Nevmian  and  P.  Tulyathan,  "Wire  Antennas 
in  the  Presence  of  a Dielectric/Ferrite  Inhomogeneity,"  IEEE 
Trans.  Antennas  and  Propagation,  Vol.  AP-26,  No.  4,  July 
1978,  pp.  587-593  (see  Appendix  A). 

Oral  Papers  E.  H.  Newman,  "Wire  Antennas  In  the  Presence  of  Ma- 
terial Bodies,"  Workshop  on  Electrically  Small  Antennas," 

Ft.  Monmouth,  N.  J.,  May  1976. 


E.  H.  Mewman,  "Wire  Antennas  in  the  Presence  of  Material  Bodies,” 
1976  IEEE/APS  International  Symposium  and  USNC/URSI  Meeting, 
Amherst,  Mass.,  October  19746. 


r 


E.  H.  Newman  and  C.  H.  Walter,  "A  Comparison  of  the  Efficiency 
of  the  Electrically  Small  Dipole  and  Loop  Antennas,"  1977  Symposium 
on  Antenna  Applications  at  Allerton  House,  Monticello,  Illinois, 
April  1977. 

Degrees  P.  Tulyathan,  M.S.  Thesis,  "Wire  Antennas  in  the  Presence 
of  a Dielectric  and/or  Ferrite  Inhomogeneity,"  The  Ohio  State 
University,  December  1976. 

P.  Mangiacotti,  M.S.  Thesis,  "The  Effects  of  the  Human  Body 
on  the  Patterns  of  the  Manpack  Transceiver  Antennas,"  The  Ohio 
Ohio  State  University,  June  1977. 

Personnel  E.  H.  Newman,  P.  Tulyathan,  C.  H.  Walter,  P.  Mangiacotti. 


Summary 


Wire  antennas  in  general,  and  electrically  small  antennas  in  par- 
ticular, often  operate  in  the  presence  of  material  bodies.  The  two 
most  common  examples  are  the  ferrite  loaded  loop  antenna  and  the  man- 
pack  whip  radiating  near  the  human  operator.  In  these  cases  the  an- 
tenna is  a thin-wire  structure  and  the  material  body  is  the  ferrite 
loading  or  the  human  operator.  This  problem  was  analyzed  via  a moment 
method  solution.  The  wire  and  material  body  were  both  represented 
by  equivalent  currents  which  were  evaluated  In  the  solution.  The  tech- 
nique and  its  associated  computer  code  are  applicable  to  thin  wire 
and  material  bodies  of  fairly  arbitrary  shape.  Further,  the  material 
bodies  can  be  lossy  and  inhomogeneous.  The  technique  was  used  to  an- 
alyze loaded  loop  antennas  and  the  effects  of  the  human  operator  on 
a manpack  antenna. 
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III.  SECOND  YEAR  - CLOSELY  SPACED  THIN-WIRES 


Journal  Publication  P.  Tulyathan  and  E.  H.  Newman,  "The  Circumfer- 
ential Variation  of  the  Axial  Component  of  Current  in  Closely 
Spaced  Thin-Wire  Antennas,"  IEEE  Trans,  on  Antennas  and  Propa- 
gation, scheduled  for  publication  January  1979  (see  Appendix 
B). 

Oral  Paper  E.  H.  Newman  and  C.  H.  Walter,  "Recent  Advances  in 
the  Method  of  Moments,"  1977  Symposium  on  Antenna  Applica- 
tions at  Allerton  House,  Monticello,  Illinois,  April  1977. 

Scientific  Personnel  E.  H.  Newman,  P.  Tulyathan,  C.  H.  Walter. 


Summary 

When  two  (or  more)  parallel  wires  come  in  close  proximity  to  each 

other,  this  causes  the  axial  current  to  be  non-uniformly  distributed 

around  the  circumference  of  the  wires.  This  so-called  proximity  effect 

2 

increases  the  ohmic  or  I R losses  of  the  wire.  In  the  design  of  the 
electrically  small  multiturn  loop  antenna  (MTL),  the  turns  are  kept 

close  together  to  minimize  the  size  of  the  antenna.  Also  the  wire 

2 

is  kept  as  thick  as  possible  to  minimize  I R losses.  Initially  as 

2 

the  wire  gets  thicker  the  I R losses  decrease  since  the  wire  current 
has  a larger  surface  area  on  which  to  spread  out.  As  wire  size  con- 
tinues to  increase  the  proximity  effect  losses  increase  to  the  point 
that  further  increasing  wire  size  actually  Increases  the  total  ohmic 
loss  and  decreases  the  antenna  efficiency.  Thus,  a knowledge  of  the 
proximity  effects  or  the  circumferential  distribtion  of  axial  current 
in  closely  spaced  wires  is  Important  to  the  design  of  the  electrically 
small  MTL. 
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The  circumferential  distribution  of  axial  current  was  studied 
using  a Fourier  Series  to  represent  the  circumferential  (or  4>)  vari- 
ation of  the  current.  In  addition  to  providing  a design  tool  applicable 
to  the  MTL,  this  study  provided  useful  information  and  insight  into 
the  thin-wire  approximation.  The  most  important  assumption  in  thin- 
wire  theory  is  that  the  axial  current  is  uniformly  distributed  around 
the  circumference  of  the  wire.  In  order  to  make  this  approximation 
it  is  required  that  the  wire  radius  a « X and  no  two  parallel  wires 
pass  within  a few  wire  diameters  of  each  other.  However,  the  results 
of  our  study  indicated  that  there  is  substantial  <{>  variation  of  the 
current  for  wire  of  radius  a = 0.001X  and  separated  by  100  wire  di- 
ameters. The  question  is  then  raised  if  the  thin-wire  assumption  and 
the  many  thin-wire  computer  codes  are  valid.  The  answer  is  that  they 
are  valid  and  useful  tools  for  predicting  such  quantities  as  far  zone 
patterns  and  input  impedance.  While  thin-wire  theory  does  not  predict 
the  proper  (J>  variation  of  the  surface  current,  it  does  predict  the 
proper  average  (over  <J>).  Since  impedance  and  far-zone  patterns  are 
primarily  dependent  upon  the  average  values  of  the  surface  current, 
thin-wire  theory  is  useful  for  computing  these  quantities. 

IV.  THIRD  YEAR  - SMALL  ANTENNA  LOCATION  SYNTHESIS 

Journal  Publication  E.  H.  Newman,  "Small  Antenna  Location  Syn- 
thesis Using  Characteristic  Modes,"  submitted  for  publication 
to  the  IEEE  Trans,  on  Antennas  and  Propagation  (see  Appendix 
C). 

Oral  Paper  C.  H.  Walter,  "Electrically  Small  Antqnna  Studies 
at  OSU,"  1978  International  Symposium  on  Antennas  and  Prop- 
agation, Sendai,  Japan,  August  1978. 

Personnel  E.  H.  Newman,  C.  H.  Walter,  P.  Tylyathan,  D.  W.  Epps 
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Summary 


The  usual  problem  with  electrically  small  antennas  is  that  they 
have  low  efficiency.  The  conventional  approaches  to  increasing  the 
efficiency  are  either  to  reduce  the  ohmic  losses  (say,  by  using  thick 
highly  conducting  wires  or  low  loss  tuning/matching  circuits)  or  to 
increase  the  radiation  (say,  by  top  loading  a dipole  or  ferrite  loading 
a loop).  Here  we  investigate  a different  approach  to  increasing  the 
radiation. 

In  actual  operation  small  antennas  are  often  mounted  on  a much 
larger  support  structure,  such  as  a jeep,  a ship,  an  airplane,  etc. 

The  support  structure  can  be  an  effective  radiator  providing  that  sub- 
stantial currents  can  be  excited  on  it  and  that  these  currents  are 
good  radiators.  The  small  antenna  is  viewed  not  as  the  primary  radiator, 
but  rather  as  a probe  to  excite  currents  on  the  support  structure. 

The  design  procedure  developed  here  employs  characteristic  modes 
to  determine  the  optimum  frequency  and  location  of  the  small  antenna 
on  the  support  structure.  In  an  example  of  a small  loop  mounted  on 
a crossed  wire  structure,  the  adiation  and  efficiency  were  increased 
by  about  a factor  of  65. 
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Wire  Antennas  in  the  Presence  of  a Dielectric/Ferrite  Inhomogeneity 

I DWARD  H.  NEWMAN,  MEMBER,  IEEE,  AND  P.  TULYATHAN 


/Itiimi  A moment  method  solution  for  treating  th in-wire 
antennas  in  the  presence  of  an  arbitrary  dielectric  and/or  ferrite  inho- 
mogeneity  is  presented.  The  wire  is  modeled  by  an  equivalent  surface 
current  density,  and  the  dielectric/ferrite  inhomogeneity  is  modeled  by 
equivalent  volume  polarization  currents.  The  conduction  currents  on 
the  wire  and  the  polarization  currents  in  the  dielectric/ferrite  inho- 
mogeneity  are  treated  as  independent  unknowns  and  determined  in 
the  moment  method  solution.  The  m»  ihod  is  applied  to  the  problem  of 
a loop  antenna  loaded  with  dielectric  or  ferrite.  Numerical  results  are 
presented,  and  are  in  good  agreement  with  measurements  and  previous 
calculations. 

I.  INTRODUCTION 

THE  INTERACTION  of  electromagnetic  fields  with  metallic 
as  well  as  dielectric  and/or  ferrite  inhomogeneities  is  a 
problem  of  practical  interest.  Numerous  techniques  have  been 

Manuscript  received  I ehruary  1 1. 1977;  revised  September  14, 1977. 
This  work  was  supported  in  pari  between  the  U.S.  Army  Research 
Oflice.  Research  Triangle  Park,  NC.  and  the  Ohio  State  University 
Research  I'outidation  under  Grant  l)AAt;29-76<i-0067. 

The  authors  arc  with  the  l.lectroScience  Laboratory,  Department  of 
Electrical  Engineering.  Ohio  State  University,  Columbus,  OH  43212. 


employed  by  various  authors  to  treat  antenna  or  scattering 
problems  with  dielectric/ferrite  inhomogeneities  [ 1 1 — C 5 ] . 
Each  of  the  above  analyses  treated  a specific  or  very  limited 
class  of  antennas  or  inhomogeneities.  Often  the  geometries 
were  two-dimensional  or  coincided  with  a separable  coordinate 
system.  It  is  desirablq  to  have  a single  method  to  analyze  a 
general  class  of  problems  where  the  antennas,  scatterers,  and 
inhomogeneities  can  be  of  arbitrary  shape,  and  the  inhomo- 
geneities can  be  dielectric  and/or  ferrite. 

The  purpose  of  this  paper  is  to  present  a moment  method 
solution  [6]  to  the  problem  of  thin-wire  antennas  and  scat- 
terers in  the  presence  of  a dielectric/ferrite  body.  The  analysis 
has  the  advantage  that  it  is  applicable  to  a wide  variety  of 
antenna,  scatterer,  and  dielectric/ferrite  geometries.  It  is  suffi- 
ciently general  to  treat  isotropic,  lossy,  and  inhomogeneous 
bodies.  The  dielectric/ferrite  inhomogeneities  are  modeled  by 
equivalent  volume  polarization  currents,  and  the  wire  is 
modeled  by  equivalent  surface  currents.  The  complex  magni- 
tudes of  these  cunents  are  determined  in  the  course  of  the 
moment  method  solution.  The  electrical  sizes  of  the  antennas, 
scatterers,  and  dielectric/ferrite  inhomogeneities  are  limited  by 
the  finite  computer  storage. 


II.  THEORY 
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Following  the  development  of  Richmond  [7]  and  Newman 
(8|,  the  moment  method  solution  for  thin-wire  antennas  and 
scatterers  in  the  presence  of  a linear  isotropic  dielectric  and/or 
femte  inhomogeneity  will  now  be  developed.  Fig.  1(a)  shows 
the  geometry  to  be  considered.  The  impressed  sources  (J,,  M,) 
are  confined  to  the  volume  Vlt  and  radiate  the  field  (E,  H)  in 
the  presence  of  two  inhomogeneities.  The  first  inhomogeneity 
is  a thin-wire  structure  in  the  volume  Vw  and  enclosed  by  the 
surface  Sw.  The  second  inhomogeneity  is  the  dielectric/ferrite 
body  in  the  volume  V2.  The  permittivity  and  permeability  of 
this  inhomogeneity  are  defined  by  the  parameters  (p , e),  which 
may  be  complex  functions  of  position.  Although  for  simplicity 
not  shown  here,  p and  < could  be  considered  as  tensor  quan- 
tities. Ibis  would  permit  the  treatment  of  anisotropic  inhomo- 
geneities.  The  ambient  medium  has  paramct-'rs  (p <>,  e0).  All 
sources  and  fields  are  considered  to  be  time  harmonic,  and  the 
e*u'  time  dependence  will  be  suppressed. 

The  first  step  m the  solution  will  be  to  replace  the  two 
inhomogeneities  by  equivalent  sources.  Employing  the  surface 
equivalence  principle  of  Schelkunoff  [9] , the  wire  is  removed 
and  the  following  surface  current  densities  are  introduced  on 
the  surface  Sw : 


(b) 

Fig.  1.  (a)  Wire  antenna  in  presence  of  dielectric  and/or  ferrite  inho- 
mogeneity. (b)  Equivalent  problem. 
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J,  =AXH  (1) 

M,  = E X n (2) 

where  n is  the  outward  directed  normal  vector  on  Su.  By 
defining  (J(,  M, ) as  in  (1)  and  (2),  the  total  field  inside  the 
wire  is  zero.  Next,  the  dielectric/ferrite  inhomogeneity  is 
replaced  by  the  ambient  medium  and  the  equivalent  volume 


polarization  currents 

J = /w(e  - «0)E 

(3) 

M =/w(p-po)H 

(4) 

TABLE  I 


field  of  t**  source  radletlM 
Source  Current  In  the  eafclent  ■edti*  (w#f«o) 
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which  are  confined  to  the  volume  V2. 

The  equivalent  problem  is  shown  in  Fig.  1(b).  Here,  in  the 
homogeneous  medium  (Po . *o)>  the  sources  (J/t  M,),  (J,,  M,), 
and  (J,  M)  radiate  the  field  (E,  H)  exterior  to.Su,  and  the  field 
(0,  0)  interior  to  Sw.  We  use  the  notation  that  the  sources 
(J<.  M,),  (J„  M,),  (J),  and  (M),  radiating  in  the  ambient 
medium  (p<>,  eo>-  produce  the  fields  (E,,  H(),  (E„  H,),  (EJ,  HJ), 
and  (EM,  Hw),  respectively.  These  relationships  are  summa- 
rized in  Table  I. 

Let  us  consider  the  wire  to  have  a circular  cross  section.  At 
each  point  on  the  wire  we  define  a local  right-handed  orthog- 
onal coordinate  system  with  unit  vectors  (A,  0,  /)  where  A is 
the  outward  normal  vector  to  Su . / is  directed  along  the  wire 

am  . a"'. 


C>  = ; / A . 


(5) 


If  the  wire  radius  a < X,  then  the  surface  current  density  on 
the  wire  structure  can  be  approximated  by  the  “thin-wire 
approximation" 


J.(f) 


?/(/)  ^ KO 

2*e  2*e 


(6) 


where  /(/)  is  the  total  current  (conduction  plus  displacement). 

The  reaction  integral  equation  (RIF.)  for  HI)  is  obtained  by 
placing  a test  source  with  current  (im,  Mm)  in  Vm,  and  noting 


that  since  a null  field  exists  in  Vw,  (Jm,  Mm)  will  have  zero 
reaction  with  the  sources  (J(,  M(),  (J,,  M,),  and  (J,  M).  Denoting 
(Em,  Hm)  as  the  fields  radiated  by  (Jm,  Mm)  in  the  ambient 
medium  (p0,  e0)  and  Z,  as  the  wire  surface  impedance  for 
exterior  excitation,  the  RIE  is 

- [L  niwr  -z.H<r)di 


-/// (J  • E"  -M  • Hm)du  ■ Vm 

(7) 

where  L represents  the  overall  wire  length  and 

Vm=  jj  j 1(1,  *E»  — M,  • Hm)de 

(8) 

*'1 

1 ft”  \ 

/ (f>E")d0 

2*  Jo 

(9) 

1 /*■  . 

(10) 

Equation  (7)  is  an  integral  equation  for  the  three  unknown 
currents  /(/),  J,  and  M.  This  equation  insures  that  these  currents 
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have  the  proper  reaction  with  a test  source  in  the  wire  surface, 
but  does  not  insure  that  the  proper  conditions  are  satisfied  in 
the  dielectnc/fernte  inhomogeneity.  In  the  dielectric/ferrite 
inhomogeneity  it  is  required  that  (E,  H)  be  related  to  (J,  M)  by 
(3)  and  (4).  Thus  we  have 


E,  + EJ  + EH  + E,  = E = 


/cu(e  - e0) 


in  V2  (11) 


H.  +H'  + HH  +Hi  = H = , in  K2.  (12) 

/GJ(  U~Mo) 


Multiplying  ( 1 1 ) by  the  vector  weighting  function  wm  ,(12) 
by  the  vector  weighting  function  wm',  and  integrating  over  the 
volume  V 2 yields 

= -/ //Ei  • wm  do  (13) 

v2 

= ~ff  f Hi'  '"m'dO.  (14) 


Equations  (7),  (13),  and  (14)  are  three  coupled  integral 
equations  which  can  be  solved  for  the  unknown  currents  /(/) 
and  (J,  M).  The  moment  method  solution  to  these  equations  is 
presented  below.  Strictly  speaking,  (7),  (13),  and  (14)  must  be 
satisfied  by  an  arbitrary  (Jm,  Mm),  wm,  and  wm'.  However,  in 
the  moment  method  solution  we  only  enforce  these  equations 
for  N distinct  ( Jm , Mm ),  M distinct  wm , and  P distinct  wm  In 
this  case  (7),  (1 3),  and  (14)  represent  N,  M,  and  P simultaneous 
linear  integral  equations,  respectively. 

The  next  step  in  the  moment  method  solution  will  be  to 
transform  these  simultaneous  linear  integral  equations  to 
simultaneous  linear  algebraic  equations.  This  is  accomplished 
by  expanding  the  unknown  currents  I,  J,  and  M as  follows: 

v 

ko  = 2 /»F»(/)  (15a) 


j=  2 /»c» 

1 

n+m+p 

M = 2 '»Qn 


In  (1  5)  the  wire  current  is  expanded  in  terms  of  the  basis  set 
F„ , and  the  electric  and  magnetic  volume  polarization  currents 
in  terms  of  the  basis  sets  G„  and  Q„,  respectively.  To  simplify 
the  following  equations  we  will  denote  (E,„,  H,„),  ( EnJ , H„0, 
and  (E„H,  H„H)  as  the  fields  radiated  by  F„,  G„,  and  Q,,, 
respectively,  in  the  homogeneous  medium  (Mo,  «<>)•  This 
notation  is  summarized  in  Table  I. 

Substituting  (13)  into  (7),  (13),  and(14),  and  changing  the 


order  of  integration  and  summation  yields 

2 [L  FmVMm  ~Z.H^)dl 

n-1  | Jo 

+ *2  Ai  l / / /Cn  * Em  dv 

n~N+ 1 l J 

v2 

N + M+P  [ f f f 

+ 2 In\+  Qn  -»mdV 

n-N+M+1  ( 

v2 

= Vm.  m = 1,2,  —,N 

Y l\  ( ( • W_  <fy| 


wm  dv 


I,  'j///E- ' w-  "I 

v2 

N*M*P  \ r r r 

+ Y / / /EnH  • WmdU 

n-N+M  + 1 

= ~fff^i’^mdv,  m 1,N+  2, -,N  + M 

(16b) 

2 ("{/// H">  ' wn ' 

+„. JL'" 

= - j f f H‘  ' w m'dv , 


m = N + M + \,N  + M + 2,  -N  + M + P.  (16c) 
Equation  ( 1 6)  can  be  written  compactly  as 

N*M*P 

2 lnzmn  = V m,  m = 1 , 2,  N + M + P (17) 


or  in  the  matrix  form  as 

Z/=V  (18) 

where  Z is  the  square  impedance  matrix,  V is  the  excitation 
voltage  column,  and  I is  the  current  column  which  contains 
the  unknown  coefficients  /„,  n = 1,  2,  •••,  jV  + Af  + P,  as 
defined  in  (13).  In  Fig.  2,  the  impedance  matrix  Z is  symbol- 
ically shown  divided  into  nine  regions,  and  the  type  of  coupling 
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Fig,  2.  Symbolic  representation  of  impedance  matrix  Z. 

accounted  for  in  each  region  is  indicated.  In  (16b)  and  (16c), 
when  m = n,  one  must  evaluate  fields  inside  a current  distri- 
bution. This  problem  is  discussed  in  the  Appendix. 
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In  order  to  obtain  numerical  results  from  the  above  formu- 
lation, it  is  necessary  to  define  specific  expansion  and  weighting 
functions.  The  choices  made  for  this  study  will  now  be  pre- 
sented. For  F„  we  choose  the  piecewise-sinusoidal  function 
used  by  Richmond  ( 7) . The  unit  magnitude  vector  volumetric 
pulse  function  is  chosen  for  the  expansion  modes  G„  and  Q„. 
The  volumetric  pulses  are  the  parallelepipeds  obtained  by 
dividing  the  dielectric/magnetic  body  into  smaller  rectangular 
volumes.  Since  the  volume  polarization  current  has  an  arbitrary 
polarization,  it  is  necessary  to  have  three  orthogonal  vector 
volumetric  pulses  occupying  the  same  parallelepiped.  Each  will 
have  a different  polarization,  Le.,  x,  y,  and  1.  There  are  no 
magnetic  test  sources  and  Mm  = 0.  The  piecewise-sinusoidal 
test  function  is  chosen  for  Jm.  The  choice  of  piecewise- 
sinusoidal  functions  for  both  the  expansion  and  test  modes 
enables  us  to  use  the  computer  program  for  thin-wire  antennas 
and  scatterers  in  a homogeneous  medium,  developed  by 
Richmond  [10],  to  calculate  the  wire/wire  region  of  the 
impedance  matrix  Z.  For  the  wm  or  wm'  we  choose  a delta 
function  which  is  located  in  the  center  of  the  corresponding 
volumetric  pulse  expansion  Gm  or  Q„.  and  with  the  same 
polarization. 

Since  the  test  modes  in  the  wire  structure  have  the  same 
current  distribution  as  the  expansion  mode,  this  may  be 
regarded  as  an  application  of  Galerkin’s  method.  The  use  of 
delta  functions  as  the  test  modes  for  the  dielectric/ferrite 
inhomogeneity  results  in  (13)  and  (14)  being  satisfied  at 
discrete  points  in  V2.  This  may  be  recognized  as  an  application 
of  the  point-matching  method. 

HI.  EXPERIMENTAL  AND  NUMERICAL  RESULTS 

In  this  section,  numerical  results  based  on  equations  derived 
in  the  preceding  section  are  presented  for  a loop  antenna 
loaded  with  a dielectric  and/or  ferrite  core.  The  calculated 
input  admittance  of  the  dielctric  loaded  loops  are  compared 
with  experimental  results.  In  the  case  of  ferrite  loaded  loops, 
the  computed  radiation  resistances  are  compared  with  the  pre- 
vious theoretical  results  by  Stewart  [1).  The  loop  geometry 
was  chosen  because  it  is  a basic  geometry  and  yet  still  required 
treating  wire  segments  which  are  parallel  or  intersect  at  an 
angle. 

The  measurements  were  made  with  a half-loop  mounted  on 
a 4-ft  square  aluminum  ground  plane.  The  half-loop  is  made  of 
copper  wire  coated  with  a thin  layer  of  tin  and  has  a total 
diameter  of  0.032  in.  The  coaxial  feed  has  an  inner  radius  of 
1/16  in  and  an  outer  radius  of  3/8  in.  The  half-loop  has  a 
dimension  of  3.02  in  by  1 .51  in.  The  dielectrics  are  rectangular 
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parallelepipeds  and  have  the  dimensions  D 1/2,  D2,  and  D3. 
The  calculations  were  made  for  the  complete  loop,  and  it  is 
the  complete  loop  admitt  ices  which  are  plotted.  The  feed 
was  modeled  by  a delta-gap  generator. 

Three  examples  will  now  be  presented  to  demonstrate  the 
ability  of  the  theory  to  treat  different  types  of  dielectric 
inhomogeneities.  Figs.  3(a)  and  (b)  compare  measured  and  cal- 
culated admittance  of  the  loop  loaded  with  Eccoflo-Hik  (er  = 
10.0,  produced  by  Emerson  and  Cuming).  The  dimensions  of 
the  block  are  D\  = 2.60  in,  D2  = 2.36  in,  and  D3  = 1.02  in. 

Next,  the  ability  of  the  theory  and  computer  program  to 
treat  problems  with  inhomogeneous  dielectrics  will  be  demon- 
strated. As  seen  in  the  insert  in  Fig.  4(a),  the  loop  is  loaded 
with  an  inhomogeneous  dielectric  of  the  same  size  as  was  con- 
sidered above.  The  dielectric  is  half-filled  with  Teflon  (er  = 
2.1)  and  half-filled  with  Eccoflo-Hik.  Figs.  4(a)  and  (b)  show 
the  measured  and  calculated  admittance  of  the  loaded  loop. 

Distilled  water  (er  = 78.05  - /3.74  at  the  resonance  fre- 
quency of  the  loop)  is  used  to  demonstrate  the  ability  to  treat 
lossy  dielectrics.  Due  to  the  large  relative  permittivity  of  water, 
the  size  of  the  inhomogeneity  has  been  reduced  to  D1  = 
1.97  in,  D2  = 1.S4  in,  and  D3  = 0.75  in.  The  measured  and 
calculated  admittance  of  the  water  loaded  loop  are  shown  in 
Figs.  5(a)  and  (b). 

One  of  the  most  common  uses  of  ferrites  is  to  increase  the 
radiation  resistance  of  loop  antennas.  Fig.  6 shows  the  radi- 
ation resistance  of  a ferrite  loaded  loop  versus  nr.  Curves  are 
shown  for  various  sizes  of  the  ferrite  cores.  The  loop  has  a 
perimeter  L = 47.24  in,  and  the  ferrite  cores  are  of  size  D 1 = 
9.45  in,  D2  = 9.45  in,  and  D3  = 2.36  in,  11.81  in,  and 
23.62  in.  The  computations  are  at  50  MHz  where  the  loop 
circumference  is  about  0.20X.  The  present  results  are  com- 
pared with  the  previous  approximate  results  of  Stewart.  These 
two  results  are  seen  to  agree  only  qualitatively.  It  is  felt  that 
most  of  the  quantitative  differences  are  due  to  approximations 
given  by  Stewart  for  the  demagnetization  factor  of  a rectan- 
gular rod.  More  recently,  Pettengill  et  at  [11]  have  presented 
different  expressions  for  the  demagnitization  factor  of  a core 
with  length  to  diameter  ratio  greater  than  two.  The  D3  = 
23.62  case  is  compared  to  a computation  using  Lcttengill’s 
demagnitization  factor,  and  the  agreement  is  good. 

IV.  SUMMARY 

The  purpose  of  this  study  has  been  to  develop  a technique 
to  analyze  thin-wire  antennas  and  scatterers  in  the  presence  of 
a dielectric  and/or  ferrite  inhomogeneity.  The  method  pre- 
sented is  a moment  method  solution  and  a modification  of  the 
piecewise-sinusoidal  reaction  formulation  for  thin-wire 
antennas  and  scatterers  in  a homogeneous  medium.  The  tech- 
nique is  sufficiently  general  to  be  applicable  to  lossy  and 
loaded  thin-wire 'antennas  and  scatterers  in  the  presence  of 
isotropic,  inhomogeneous,  and  lossy  dielectric/ferrite 
inhomogeneities. 

Numerical  results  have  been  presented  and  compared  with 
measurements  and  previous  calculations.  These  results  show 
the  ability  of  the  theory  and  computer  program  to  analyze 
antennas  in  the  presence  of  dielectrict  of  varying  sizes  and 
permittivities,  inhomogeneous  dielectrics,  lossy  dielectrics,  and 
ferrites.  The  method  is,  however,  limited  by  the  finite  com- 
puter core  to  electrically  small  inhomoneneiiies. 

Using  the  technique  presented  here  for  treating  antennas 
and  scatterers  in  the  presence  of  a dielectric /ferrite  inhomo- 
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Fig.  3.  (a)  Comparison  of  calculated  and  measured  conductance  for 

dielectric  («r  * 10.0)  loaded  loop,  (b)  Comparison  of  calculated  and 
measured  susce,  ance  for  dielectric  <tr  » 10.0)  loaded  loop. 

Fig.  4.  (a)  Comparison  of  calculated  and  measured  conductance  for 
dielectric  (inhomogeneous)  loaded  loop,  (b)  Comparison  of  cal- 
culated and  measured  susceptance  for  dielectric  (inhomogeneous) 
loaded  loop. 
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Fig.  6 Comparison  of  present  theory  with  previous  calculation  for 
the  radiation  resistance  of  ferrite  loaded  loops. 


geneity,  one  can  analyze  problems  such  as  ferrite  loaded  loops, 
manpack  transceiver  antennas,  scattering  from  rain  or  blood 
cells,  effects  of  microwave  radiation  on  biological  tissues,  and 
antennas  covered  by  dielectric  radomes. 


APPENDIX 


Consider  the  problem  of  finding  the  electric  field  at  the 
center  of  a column  of  current  with  uniform  current  density 
J = zJ,  occupying  a circular  cylindrical  volume.  In  numerically 
evaluating  this  field,  one  must  exclude  a small  region  about  the 
field  point  Unfortunately,  the  value  of  the  resultant  integral  is 
dependent  upon  the  shape  of  the  volume  excluded.  If  one 
excludes  a circular  cylinder  of  radius  r)  and  height  2 h sur- 
rounding the  singularity,  Van  Bladel  [12]  found  that  the 
result  will  he  in  error  by 
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Note  that  if  one  excludes  a long  thin  cylinder  where  t}<h,6 
goes  to  zero  and  the  correct  result  is  obtained.  This  idea  is 
used  in  determining  the  electric  field  inside  a volumetric 
electric  current  expansion  mode,  and  also  for  the  dual  problem 
of  determining  the  magnetic  field  inside  a volumetric  magnetic 
current  expansion  mode. 


100*  ««IM1U  (*y\  ) 

(b) 

Fig  J (a)  Comparison  of  calculated  and  measured  conductance  for 
dklectnc  (distilled  water)  loaded  loop,  fb)  Comparison  of  cal- 
culated and  measured  suaceptance  for  dielectric  (distilled  water) 
loaded  loop. 
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THE  CIRCUMFERENTIAL  VARIATION  OF  THE  AXIAL  COMPONENT  OF 
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ABSTRACT 

Closely  spaced  thin-wire  antennas  are  analyzed  by  the  moment  method 
using  the  piecewise-sinusoidal  function  to  describe  the  current  variation 
along  the  length  of  the  wire  and  a Fourier  series  for  the  circumferential 
variation.  Data  are  presented  for  the  circumferential  variation  of  the 
surface  current  density  of  parallel  dipoles  obtained  by  the  "thin-wire" 
theory,  wire-grid  model,  and  the  present  formulation.  It  is  found  that 
there  is  substantial  circumferential  variation,  even  for  dipoles  spaced 
more  than  several  wire  diameters. 
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DAAG29-76-G-0067  between  the  Department  of  the  Army,  U.  S.  Army  Research 
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I.  INTRODUCTION 

A fundamental  approximation  made  in  the  "thin-wire"  theory  is  that  the 
surface  current  density  is  uniform  around  the  circumference  of  the  wires, 
in  order  to  make  this  approximation,  two  conditions  are  imposed  on  the  wires. 
First,  the  wire  radius  must  be  much  less  than  a free  space  wavelength,  and 
second,  no  two  wires  pass  within  a few  wire  diameters  of  each  other. 

This  paper  considers  the  problem  of  closely  spaced  three-dimensional 
electrically  thin  wires.  Smith  [1]  and  Olaofe  [2]  had  considered  the  problem 
of  closely  spaced  infinite  cylinders,  and  Smith  included  a modification  for 
finite  wires.  The  solution  presented  here  is  a modification  of  the  piecewise- 
sinusoidal  reaction  formulation  for  "thin  wire"  structures  [3].  A Fourier 
series  is  used  to  represent  the  circumferential  variation  of  the  axial  com- 
ponent of  the  current.  Numerical  data  illustrate  the  effects  of  wire  radius 
and  separation  on  the  circumferential  variation  of  the  current.  The  data 
are  compared  with  the  results  of  "thin-wire"  theory  and  also  with  a wire-grid 
model  of  the  closely  spaced  wires. 

II.  THEORY 

The  reaction  concept  was  introduced  by  Rumsey  [4],  and  applied  by 
Richmond  [3]  to  derive  the  reaction  integral  equation  (RIE)  for  radiation 
and  scattering  from  thin-wires  and  conducting  surfaces. 

Consider  the  problem  of  scattering  by  a wire  structure  enclosed  by 
the  surface  S,  and  with  the  impressed  electric  and  magnetic  currents 
confined  in  the  volume  V.  Let  the  ambient  medium  be  free  space.  The 
e^wt  time  dependence  of  all  sources  and  fields  will  be  suppressed.  In  the 
presence  of  the  wire,  generate  the  field  (E..H.).  From  the  surface- 

equivalence  theorem  of  Schelkunoff  [5],  the  wire  can  be  replaced  by  free  space 

if  the  following  surface-current  densities  are  Introduced  on  the  surface  S 
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^ = n x H 


(1) 


l 


Hj’Exn  (2) 

where  n is  the  outward  directed  unit  vector  on  S.  By  defining  (Jg.Mg)  as  in 
Equations  (1)  and  (2),  the  total  field  inside  the  wire  is  zero.  The  RIE  for 
the  unknown  currents  is  [3] 

j]  (4*im-Ms'Hm)ds  + jjj  (Ji*£m-M.*Hm)dv  = 0 (3) 

s v 

where  (Em,H_m)  is  the  field  of  test  sources  located  interior  to  the 

surface  S and  radiating  in  the  homogeneous  medium.  Below  we  will  consider 
perfectly-conducting  wires,  and  thus  * 0. 

Let  each  segment  of  the  wire  structure  have  a circular  cross  section 
and  at  each  point  on  the  surface  define  a right-handed  orthogonal  coordinate 
system  with  unit  vectors  (n,4,i)  where  n is  the  outward  normal  vector,  £ is 
directed  along  the  wire  axis  and 

A A A 

♦ a l x n.  (4) 

At  this  point.  Equation  (3)  is  usually  simplified  by  making  the 
following  thin-wire  approximations: 

1.  wire  radius  a « A 

2.  neglect  integrations  over  the  end  surfaces  of  the  wire 

3.  the  ♦ component  of  is  zero 

A 

4.  the  a component  of  is  Independent  of  ♦. 

We  will  use  all  of  the  above  approximations  except  for  #4.  The  surface 
current  density  on  the  wire  can  then  be  written  as 

«>,(*.♦)  - * JS(M) 


F 

fl 
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(5) 


Inserting  Equation  (5)  into  Equation  (3)  yields 


L 2it 

- j j 0S(M)  (*-£m)  ad*  d*  - jjj  CJi*£m-Mi*!im)dv  , (6) 

0 0 V 


where  L denotes  the  overall  wire  length. 

Equation  (6)  will  be  solved  via  the  method  of  moments  [6].  To  do 
this  Js (*■  ,4> ) is  expanded  in  terms  of  a finite  series  as  follows: 

N 

JS(M)  = l Jn  Gn(M)  (7) 


where  the  J are  unknown  coefficients  and  the  G„(a,4>)  are  a known  basis 
n n' 

set.  Enforcing  Equation  (6)  for  N distinct  test  sources  yields  the  following 
system  of  simultaneous  linear  equations 


N 

l Jn  Zmn  = Vm 
n=i  n mn  m 


m = 1 ,2,*«  *N 


where 


Zmn  = -|j  G„(A^)  (£- f)  & (9a) 

n 

Vm  - jjj  (Ji-£m-Mi-Hm)dv  (9b) 

v 


and  where  the  integration  in  Equation  (9a)  is  over  the  surface  of  the  n-th 
expansion  mode.  The  accuracy  and  computational  efficiency  of  the  solution 
are  dependent  upon  the  choices  for  the  expansion  and  testing  functions. 

In  defining  expansion  functions  we  choose  the  plecewise-sinusoldal 
function  [3]  to  describe  the  a variation  and  a Fourier  series  for  the  $ 
variation.  The  expansion  Is 
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n 

w 


T 


P K 

J«(M)  » I FU)  + l < 


P-1 


P P 


k=1 


♦ J2kp»p  y*> s,n  k*] 


it. 


)P+P  Fp(4>  cos  k* 


(10) 


Thus,  the  Gn(i,$)  are  of  the  form  Fn(ft),  FnC«- )cos<P , Fn(ft)sin<j>,  Fn(*' )cos2<|> , 
Fn(l)sin2$,  •••  Fn(4)cosK(fi  and  Fn(ft)sinK<)>  where  the  Fn(ft)  are  the 
piecewise-sinusoidal  V-dipoles.  If  P V-dipoles  are  used  to  describe  the  ft 
variation  and  K terms  in  the  Fourier  series  are  retained  .then  the  total 
number  of  unknowns  will  be 

N = P(2K+1 ) . (11) 


The  testing  functions  were  chosen  Identical  to  the  expansion  modes. 

A 

Thus,  Mp,  = 0 and  ^ = ftGm«  This  is  an  application  of  Galerkin's  method 
and  the  impedance  matrix  will  be  symmetric.  Our  initial  choice  for  the 
test  modes  was  filamentary  sources  on  the  wire  axis.  This  choice  led  to 
serious  problems  of  relative  convergence  [7]  and  was  abandoned  for  the  true 
Galerkin  solution. 

The  above  choice  of  expansion  and  testing  functions  has  several 
advantages.  First,  the  expansion  modes  are  placed  in  an  overlapping  array 
on  the  wire  so  that  continuity  of  current  is  enforced.  Evaluation  of  the 
elements  in  the  Impedance  matrix.  Equation  (9),  requires  a quadruple 
integration,  i.e.,  two  to  find  and  two  to  Integrate  over  the  surface  of 
the  n-th  expansion  mode.  However,  since  we  employ  the  piecewise-sinusoidal 
functions,  only  two  of  these  Integrations  are  done  numerically.  Advantage 
can  also  be  taken  of  the  orthogonality  of  the  Fourier  modes  In  evaluating 
terms  where  the  expansion  and  test  modes  overlap. 
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III.  NUMERICAL  RESULTS 

In  this  section  numerical  results  will  be  presented  for  the  surface 

current  density  of  parallel  center-fed  dipoles,  as  illustrated  in  Figure  1. 

Results  obtained  by  the  method  of  the  previous  section  will  be  compared 

to  those  of  "thin-wire"  theory  [8],  and  also  to  the  results  of  a wire-grid 

model  of  the  parallel  dipoles  [9].  In  the  wire-grid  model  each  thin-wire 

dipole  is  modeled  by  about  sixteen  extremely  thin  wires  uniformly  spaced 

around  its  circumference.  The  dipoles  are  center-fed  by  in  phase  unit 

voltage  delta-gap  generators  at  30  MHz.  Although  varies  along  the  t 

coordinate,  only  J around  the  center  of  the  wire  will  be  presented.  The 
«> 

following  calculations  demonstrate  the  ability  of  the  preceding  formulation 
and  computer  programs  to  account  for  the  <j>  variation  of  for  closely  spaced 
thin  and  moderately  thick  parallel  dipcles.  The  surface  current  densities 
are  normalized  with  respect  to  an  isolated  dipole.  In  all  the  computations 
to  follow,  3 modes  were  used  to  describe  the  longitudinal  distribution,  on 
each  a/2  dipole. 

The  normalized  magnitude  and  phase  of  the  surface  current  density  of 

two  parallel  a/2  dipoles  for  various  wire  spacings,  c/a,  and  wire  radii  of 
-3  -2 

10  A and  10  a are  shown  in  Figures  2 and  3,  respectively.  Referring  to 
Equation  (10),  the  coefficients  of  the  expansion  modes  for  the  above  cases, 
and  for  a few  cases  of  very  large  c/a,  are  presented  in  Tables  1-A  and  1-B. 
Note  that  since  is  an  even  function  in  <j>  all  sink<j>  terms  vanish.  Following 
the  notation  of  Equation  (10),  K will  denote  the  number  of  <(>  dependent  terms 
in  the  surface  current  expansion.  For  example,  K=0  implies  that  only  the 
constant  ("thin-wire"  approximation)  term  is  used,  K=1  implies  that  the  con- 
stant, cos$  and  sin4>  terms  are  used,  etc.  In  Figures  2 and  3 note  that  the 
results  of  the  previous  section  are  in  excellent  agreement  with  the  wire-grid 


I 

I 


model  of  the  parallel  dipoles.  In  Tables  1-A  and  1-B  note  that  the 
coefficients  of  the  constant  term  for  the  K=0  case  are  very  close  to  those 

O 

of  the  K=1  and  K=2  cases,  except  in  the  extreme  case  where  a = 10"  X and 
c/a  =1.1. 

Figure  4 and  Figure  5 show  the  normalized  magnitude  and  phase  of  the 

surface  current  distribution  for  three  parallel  x/2  dipoles  for  various 

3 2 

wire  spacings  and  wire  radii  of  10  A and  10  A,  respectively.  The  coef- 
ficients of  the  expansion  modes  are  tabulated  in  Tables  2-A  and  2-B  for  the 
above  cases.  In  the  cases  where  the  wires  are  very  close  together  (c/a=l.l 
and  c/a=2.0),  the  coefficients  predicted  by  the  case  K=0  can  be  in  error 
by  an  order  of  magnitude. 


i 


f] 

n 
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TABLE  1-A.  Two  Parallel  A/2  Dipoles  Each  of  Radius  10"^A 


Coefficients  of  the  Expansion  Modes  for  J at  the 
Center  of  the  Dipole  #1  (mA/neter) 

Note  that  sink<f>  terms  vanish 


c/a 

expansion  mode 
used  (K) 

Constant 

cos$ 

cos2$ 

0 

86. 50/ -26. 6 

- 

• 

1.1 

1 

86.33/-26.5 

166.29/149.4 

- 

2 

86.33/-26.5 

164.65/149.4 

34.12/152.0 

0 

86.06/-26.1 

- 

2.0 

1 

86.00/-25.9 

168.12/148.2 

- 

2 

86.00/-25.9 

167.63/148.2 

30.36/149.8 

0 

! 

88.23/-20.2 

- 

- 

20.0 

1 

88.24/-20.1 

166.23/143.8 

- 

2 

88.24/-20.1 

166.23/143.8 

4.11/150.6 

0 

97.59/-13.6 

- 

- 

50.0 

1 

97.64/-13.6 

164.11/139.0 

- 

2 

97.64/-13.6 

164.11/139.0 

1.86/156.6 

0 

161.38/01 .8 

1 ,000 

1 

161.06/01.8 

46.23/-136.8 

- 

2 

161.06/01.8 

46.23/-136..P 

0.26/-62.4 

0 

1 172.97/-27.7 

- 

- 

10,000 

1 

172.64/-27.7 

5.02/-127.2 

- 

2 

172.64/-27.7 

5.02/-127.2 

0.03/-1 12.2 

0 

1 74.03^-27.7 

- 

- 

00,000 

1 

173.73/-27.2 

0.51/-170.5 

- 

2 

173.73/-27.2 

0.51/-1 70.5 

0.0 

used  (K) 


Constant 


7.96/-15.7 


8.19/-13.7 


8.19/-13.6 


8.03/-1 1 .2 

8.327-9.9 

8.327-9.9 


18.38/7.0 

18.35/6.7 

18.35/6.7 


13.49/-26.9 


13.56/-26.4 


13.56/-26.4 


16.06/-20.3 


16.21/-19.7 


16.21 /-I 9.7 


16.18/-19.8 


16. 33/-1 9.2 


16. 33/— 1 9.2 


20.64/149.0 
20.49/148.8 


21.47/145.5 


21.40/145.5 


24.51/83.1 


24.51/83.1 


10.90/-146.6 
10.90/-146.6 


0.66/-123.0 


0.66/-123.0 


5.10/157.4 


4.89/153.7 


1.93/133.6 


0.82/-76.1 


0.05/-36.5 


0.07/-166.3 
0.07/-166.3  0.0 


TABLE  2-A.  Continued 


wire  #1  <20.0 


wire  #2  <20.0 


wire  #1  <50.0 


wire  #2  <50.0 


expansion  mode 
used  (K) 


Constant  cos$ 

91.89/-21.2 

90.07/-20.9  170.14/143.7 

90.75/-20.9  170.14/143.7 

90.75/-20.9  170.14/143.7 

24.35/74.1 
23.69/68.7 
23.69/68.7  0.0 

23.69/68.7  0.0 


I03.08/-13.9 
I03.67/-13.7 
102 .67/— 1 3.7 
102.67/-13.7 

47.96/53.9 

47.96/52.7 

47.96/52.7 

47.96/53.7 


cos2$ 


cos3<j> 


2.51/160.0 
2.51/160.0  0.0 


8.36/150.2 
8.36/150.2  0.0 


169.30/138.4 

169.30/138.4 

169.30/138.4 


1.51/173.1 
1.51/173.1  0.0 


3.86/156.5 
3.86/156.5  0.0 


* ■*)'  . < »'?■  £ ■ 

. V . 7 r-'-  7 


xpansion  mode 


wire  #1(20.0 


wire  #2<20.0 


wire  #1 <50.0 


wire  #2(50.0 


■I 


TABLE  2-B.  Continued 


Constant 


17.63/2.0 

17.76/2.5 

17.76/2.5 

17.76/2.5 


24.16/21.7 


23.83/22.1 


23.83/22.1 


23.83/22.1 


cos# 

cos  2# 

20.69/78.2 

- 

20.69/78.2 

1.91/124.9 

20.69/78.2 

1.91/124.9 

11.18/-34.8 


11 .14/-34.6 


11 .12/-34.6 


11.12/-34.6 


cos  3# 


129.0  I 0.0 


12. 40/ -28. 9 

- 

- 

12.50/-28.4 

14 . 54/— 1 49 . 2 

• 

- 

12.50/-28.4 

14.54/-149.2 

1.09/-75.6 

• 

12.50/-28.4 

14. 54/- 1 49 . 2 

1.09/-75.6 

0.03/-8.6 

1.55/-78.8 
1. 5S/-78. 8 


The  results  shown  in  Figures  2-5  and  in  Tables  1 and  2 suggest  opposite 


‘ 


conclusions  as  to  the  accuracy  of  the  usual  thin-wire  approximation  that  is 
independent  of  <j>.  First,  Figures  2-5  clearly  show  that  is  strongly  $ 
dependent,  even  for  c/a  as  large  as  50  and  a as  small  as  A/1 000 . Thus,  one  is 
led  to  the  conclusion  that  the  thin-wire  approximation  is  not  valid.  Next, 
Tables  1 and  2 show  that,  except  for  extremely  small  c/a  where  the  thin-wire 
approximation  is  known  to  fail,  the  <t>  independent  term  predicted  from  "thin- 
wire"  theory  (K=0  case)  is  essentially  identical  to  the  <t»  independent  term 
from  the  K=1  or  K=2  cases.  Since  for  wires  whose  radius  a <<*  it  is  the  $ 
independent  term  which  dominates  the  far-zone  radiated  fields,  "thin-wire" 
theory  will  predict  the  proper  far-zone  fields.  Also,  if  the  device  used  to 
measure  input  impedance  is  sensitive  to  the  average  value  of  voltage/J$  around 
the  circumference  of  the  wire,  then  it  will  also  depend  only  on  the  $ inde- 
pendent modes.  Thus,  one  can  be  led  to  the  conclusion  that  the  "thin-wire" 
approximation  is  in  fact  valid.  This  discrepency  can  be  resolved  by  stating 
that  the  "thin-wire"  approximation  is  not  valid  for  predicting  the  circumfer- 
ential distribution  of  current  on  a general  wire  structure.  However,  it  is 
useful  for  computing  such  quantities  as  impedance  and  far-zone  fields. 

At  this  point  it  is  reasonable  to  ask  the  question  "why  go  to  all  the 
time  and  trouble  to  determine  the  <t>  dependence  of  ^ since  it  ordinarily 
does  not  affect  either  the  input  Impedance  or  far-zone  fields"?  The  answer 
is  two-fold.  First,  "thin-wire"  computer  programs  are  In  such  widespread 
use  that  any  information  relative  to  the  basic  approximations  made  in  these 
codes  is  of  value.  Secondly,  there  are  quantities  which  the  <t>  dependence 
of  does  affect.  For  example,  the  conductor  loss  resistance  Is  a strong 
function  of  the  circumferential  distribution  of  current  on  the  wire  [1].  A 
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knowledge  as  to  how  this  loss  resistance  varies  with  wire  radius  and  wire 
spacing  is  important  in  determining  the  efficiency  of  such  antennas  as  the 
electrically  small  multiturn  loop  antennas.  The  41  dependence  of  can 
also  be  important  to  the  transient  behavior  [10], 

IV.  CONCLUSION 

This  paper  has  considered  the  problem  of  determining  the  circumferential 
current  distribution  on  closely  spaced  electrically  thin-wires.  The  solution 
is  a modification  of  the  piecewise-sinusoidal  reaction  formulation  for  thin- 
wire  structures.  A Fourier  series  is  used  to  represent  the  <f>  dependence  of 
the  wire  surface  current  density.  Numerical  data  presented  illustrate  that 
there  can  be  a considerable  circumferential  variation  even  when  the  wire  to 
wire  separation  exceeds  several  wire  diameters.  This  can  be  important  in 
determining  antenna  efficiency;  however,  if  these  wires  were  to  be  analyzed 
by  "thin-wire"  theory  (i.e.,  neglecting  the  circumferential  variation  of  the 
surface  current  density)  the  correct  impedance  and  far-zone  fields  are 
usually  obtained. 
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APPENDIX  C 

SMALL  ANTENNA  LOCATION  SYNTHESIS 
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ABSTRACT 

It  is  shown  that  the  efficiency  of  a small  antenna  can  be  substantially 
increased  by  properly  locating  it  on  its  support  structure.  Characteristic 
modes  are  used  to  determine  the  optimum  location  and  frequency. 
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I.  DISCUSSION 

An  important  problem  in  antenna  design  is  the  development  of  rela- 
tively efficient  electrically  small  radiating  elements.  The  antenna  ef- 
ficiency can  be  defined  as 


E 


R +R« 
r * 


(1) 


where  R is  the  radiation  resistance  and  R.  is  the  loss  resistance.  Small 
r 8, 

antennas  are  generally  very  inefficient  because  their  radiation  resistance 

2 

is  low.  The  radiation  resistance  of  a small  dipole  varies  as  U/x)  , while 
that  of  a small  loop  varies  as  (fc/A)^  where  t is  the  maximum  extent  of 
the  antenna  and  X is  the  wavelength.  There  are  two  conventional  approaches 
to  improving  the  efficiency  of  small  antennas.  The  first  is  to  reduce 
the  loss  resistance,  say  by  using  thick  highly  conducting  wires  to  con- 
struct the  antenna,  or  by  using  low  loss  components  in  the  tuning  or  matching 
networks.  The  second  is  to  increase  the  radiation  resistance,  say  by  top 
loading  a short  dipole  or  ferrite  loading  a small  loop.  Below  is  described 
a third  technique  for  increasing  the  radiation  resistance. 

Small  antennas  are  usually  designed  on  test  beds  resembling  either 
free  space  or  an  infinite  ground  plane.  However,  in  use  they  are  often 
mounted  on  support  structure  such  as  a ship,  a tank,  a man,  or  an  airplane. 
The  basic  idea  here  is  to  think  of  the  small  antenna  not  as  the  primary 
radiator,  but  rather  as  a probe  to  excite  currents  on  the  support  struc- 


ture. Since  the  support  structure  often  is  not  electrically  small  it  can 
be  an  effective  radiator.  Thus,  the  radiation  resistance  and  efficiency 
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One  problem  is  where  to  locate  the  small  antenna.  Some  locations 
may  result  in  substantial  improvements  in  efficiency,  while  others  may 
result  in  little  or  no  improvement.  A second  and  related  problem  is  the 
selection  of  an  operating  frequency  which  optimizes  the  efficiency.  The 
optimum  location  at  one  frequency  may  not  be  the  optimum  at  another  fre- 
quency. The  antenna  designer  may  not  have  complete  freedom  in  selecting 
location  and  frequency.  However,  knowledge  as  to  how  the  efficiency  varies 
with  frequency  and  location  can  help  in  the  selection  of  the  best  allowable 
optimum. 

In  order  to  select  the  optimum  location  and  frequency  one  could  com- 
pute and  plot  a family  of  curves  of  efficiency  versus  position  versus  fre- 
quency. Presented  here  is  a technique  which  should  involve  considerably 
less  computation,  and  also  give  more  physical  insight.  The  technique  in- 
volves the  computation  of  the  characteristic  modes  of  the  support  structure. 
The  theory  of  characteristic  modes  has  been  presented  by  Garbacz  [Y]  and 
by  Harrington  and  Mautz  QQ  . Numerically  efficient  techniques,  using  the 
method  of  moments,  for  computing  characteristic  modes  have  also  been  pre- 
sented by  Harrington  and  Mautz  QQ,  and  are  applied  in  this  work. 

The  characteristic  modes  are  real  currents  on  the  surface  of  a con- 
ducting body.  Denoting  as  a characteristic  mode,  the  choice  as 

a basis  set  for  the  current  diagonalizes  the  impedance  matrix  or  operator 
of  the  conducting  body.  The  characteristic  modes  have  orthogonality  of 

the  radiated  fields.  Associated  with  each  characteristic  mode  J is  a 

“11 

real  characteristic  value  or  eigenvalue  A . 

The  eigenvalues  are  important  because  they  tell  how  well  a particular 
mode  radiates.  Those  modes  with  small  |AJ  are  good  or  effective  radiators, 
while  those  with  large |Xr| are  poor  or  ineffective  radiators.  In  order 
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to  substantially  improve  the  efficiency  of  a small  antenna  mounted  on  a 

support  structure,  it  is  necessary  to  excite  modes  which  are  effective 

radiators,  i.e.,  those  with  small  |Xn|  . Associated  with  each  eigenvalue 

is  a characteristic  angle  defined  by 

a = 180°  - tan  *(A  ).  (2) 

n n 

Modes  with  characteristic  angles  near  180°  are  effective  radiators,  while 
those  with  characteristic  angles  near  9C°  or  270°  are  ineffective  radiators. 

A source  or  probe,  with  impressed  field  E*-,  excites  the  nC^  charac- 
teristic mode  with  strength 


'» ■ \\  ^ 


where  the  integral  extends  over  the  surface  of  the  body.  Equation  (3) 
shows  that  in  order  to  excite  ^ as  strongly  as  possible,  the  probe  should 
be  placed  at  or  near  the  maximum  of  J^.  Further,  the  probe  should  be  oriented 
so  that  E1  and  ^ are  parallel.  Thus,  to  improve  the  efficiency  of  a small 
antenna  we  wish  to  locate  it  on  its  support  structure  where  a characteristic 
mode,  with  characteristic  angle  near  180°,  is  maximum.  The  design  example 
below  will  illustrate  this  procedure  for  a small  loop  on  a crossed  wire. 


II.  DESIGN  EXAMPLE 

A design  example  will  now  be  presented  to  illustrate  the  use  of  charac- 
teristic modes  to  select  the  operating  frequency  and  location  for  a small 
loop  on  a crossed  wire.  The  crossed  wire  is  shown  in  the  insert  in  Figure 
1,  and  could  represent  a crude  model  for  an  airplane  shape.  As  a design 
restriction  we  will  assume  that  the  loop  must  be  located  on  the  longer 
vertical  wire  of  length  L. 
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The  design  is  begun  by  computing  the  characteristic  modes  and  charac- 
teristic angles  of  the  crossed  wire  versus  frequency.  Figure  1 shows  a 
plot  of  the  characteristic  angles  of  the  first  few  modes  versus  L/X . The 
optimum  frequencies  are  those  where  the  various  modes  are  resonant,  i.e., 
a = 180°.  Figure  2 shows  a plot  of  the  characteristic  modes  at  L/X  = 0.75 
where  mode  C is  nearly  resonant.  In  this  figure  the  solid  line  represents 
the  current  on  the  vertical  wire  of  length  L,  while  the  dotted  line  represents 
the  current  on  the  horizontal  wire  of  length  2L/3.  The  arrows  indicate 
the  direction  of  current  flow.  The  characteristic  angles  of  the  modes 
are  also  shown.  Note  that  modes  A and  D have  angles  near  90°  and  27 0°, 
respectively.  Thus,  they  are  poor  or  ineffective  radiators  and  we  have 
no  interest  in  exciting  them.  Mode  B,  with  a = 148°,  is  a reasonably  good 
radiator,  but  it  has  zero  current  on  the  vertical  wire  where  the  small 
antenna  must  be  located.  Thus,  we  will  not  be  able  to  excite  mode  B with 
substantial  strength.  Finally,  mode  C is  an  excellent  radiator  with 
a = 178°.  It  has  maximum  current  just  above  the  wire  junction,  and  thus 
this  is  the  optimum  location  of  the  small  loop. 

Figure  3 shows  a plot  of  the  radiation  resistance  of  a small  loop, 

normalized  to  its  value  in  free  space,  versus  the  loop  position  on  the 

crossed  wire  at  L/X  ■ 0.75.  The  ratio  R /R  , where  R is  the  radiation 

r ro  ro 

resistance  of  the  loop  in  free  space  and  Rf  is  the  radiation  resistance 
of  the  loop  on  the  crossed  wire,  is  equal  to  the  increase  in  efficiency 
due  to  the  crossed  wire,  provided  the  efficiency  is  low.  Note,  that  as 
predicted  above,  the  radiation  resistance  increases  substantially  (by  a 
factor  of  about  65)  when  it  is  located  just  above  the  junction.  On  the 
other  hand,  if  a location  below  the  junction  (s<1.0)  is  chosen  then  the 
increase  is  less  than  about  a factor  of  two. 
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III.  SUMMARY 

It  has  been  demonstrated  that  the  radiation  resistance  and  efficiency 
of  a small  antenna  can  be  substantially  increased  by  properly  locating 
it  on  its  support  structure.  Characteristic  modes  are  shown  to  be  a useful 
tool  for  determining  the  optimum  location  and  operating  frequency.  For 
two  reasons,  the  technique  is  best  applicable  to  support  structures  which 
are  not  electrically  large.  First,  the  size  of  the  moment  method  impedance 
matrix,  necessary  for  the  computation  of  the  characteriatic  modes,  can 
become  unpractically  large.  Second,  as  the  electrical  size  of  the  support 
structure  increases,  more  and  more  characteristic  modes  are  effective 
radiators.  In  this  case,  the  exact  antenna  location  and  operating  fre- 
quency is  likely  to  be  far  less  critical. 
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